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Localization properties of random-mass Dirac fermions for a realization of mass disorder, com- 
monly referred to as Cho-Fisher model, is studied on the D-class chiral network. We show that 
a simple RG description captures accurately three phases: thermal metal and two insulators with 
quantized Hall conductances, as well as transitions between them (including critical exponents). 
We find that, with no randomness in phases on the links, transmission via the RG block exhibits a 
sizable portion of perfect resonances. Derealization occurs by proliferation of these resonances to 
larger scales. Evolution of the thermal conductance distribution towards metallic fixed point is syn- 
chronized with evolution of signs of transmission coefficients, so that derealization is accompanied 
with sign percolation. 
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Introduction. As it was pointed out in Ref. [1], there ex- 
ists a profound connection between 2D electron motion in 
random potential, V(x, y), in a quantizing magnetic field 
and the motion of 2D Dirac fermions with random mass, 
M(x, i/), in a zero field. This connection can be traced to 
the analogy of the drift of electron along closed equipo- 
tentials, V(x, y) = E, and chiral motion of Dirac fermion 
along contours, M(x,y) = 0. Due to this similarity, 
conventional quantum Hall transition, which takes place 
as equipotentials merge, has its counterpart for Dirac 
fermions. However, critical behaviors for these two tran- 
sitions are different due to the different underlying sym- 
metries. In simple terms, the difference comes from the 
fact that electron drifting along equipotential acquires a 
random phase, while Dirac fermion moving along a line 
M = does not. According to general classification 
derealization transition of Dirac fermions with mass ran- 
domness belongs to the D-class time-reversal and 
spin-rotational symmetries are broken. It was realized 
that derealization transition in the D-class is related to 
criticality in the random-bond Ising model lHZL ther- 
mal quantum Hall effect in superconductors U [9, 14-17) 
with either spin-orbit scattering or with special type of 
pairing, and lately, to in-gap transport in graphene with 
broken sublattice symmetry To|-13|, 

Not only the critical exponent of localization radius 
for the D-class is different from the quantum Hall critical 
exponent, but the very picture of derealization in the D- 
class is much less intuitive than in the quantum Hall tran- 
sition. This is because transformation of the wave func- 
tion corresponding to the drift along V(x, y) = E upon 
merging of two equipotentials is different from transfor- 
mation of chiral Dirac wave functions upon merging of 
two contours M(x, y) = 0. In both cases, the reflection 
and transmission coefficients, r and t, at the point of 
merging can be chosen real with r 2 + t 2 = 1, ensuring 
the current conservation. In the quantum Hall case, the 
signs of r and t can be chosen arbitrary, since change of 
the sign can be absorbed into the random drift-phase. 
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FIG. 1: (Color online) CF model. The boundary on the plane 
gap-disorder ((t 2 ), W) (solid line) separates metallic and insu- 
lating fixed points of the RG transformation Eq. Q. Dashed 
line-the boundary between insulators with different values of 
quantized thermal Hall conductance. Purple triangle marks a 
tricritical point at Wtr ~ 0.06. Upper and lower insets show 
the evolution of the portion, k, of negative values of reflection 
coefficients, n, with the sample size calculated for (t 2 ) = 0.2 
in metallic and insulating phases, respectively. 



By contrast, there is no such freedom in the D-class. In- 
deed, the phase accumulated by Dirac edge state along 
a contour M(x, y) = is always equal to 7r, which is the 
consequence of the pseudospin structure of the eigenfunc- 
tion [1]. If scattering matrices 



Si = 



-ti 



(1) 



are the same for all points of contact i (without any 
randomness in signs), it is obvious that derealization 
will take place at r = t = 1/V% which corresponds to 
(M(x,y)} = 0. Then the critical exponent is v = 1, 
which is a consequence of in-gap tunneling of the Dirac 
fermion. 

An amazing property of the D-class is that disorder in 
signs of and U has a drastic effect on derealization 
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FIG. 2: (Color online) Left: Chiral network for the CF model. 
Sign arrangement of t and r at the nodes ensures tt flux per 
plaquette [14] . Red lines single out a group of five nodes con- 
stituting an RG block. Right: Illustration of the RG trans- 
formation Eqs. (|4j), (|9j) for elements of renormalized node 
scattering matrix on the red sublattice. Truncation proce- 
dure (dashed lines on the original lattice) enforces, via sign 
structure in the nodes, the factors 1 and —1 in the corners. 
These factors are taken into account in Eqs.j4}, ([9]). 



transition [Hj]. Moreover, for a given degree of the net 
sign disorder, a particular way in which it is introduced 
ja 0,0,111 Hg can turn the system either into metal (M) 
or into insulator (I). In particular, if disorder is intro- 
duced by randomly changing signs of columns in Eq. (pQ) 
with probability, W, then arbitrarily small W leads to de- 
localization at any average mass, {M(x,y)) oc (£?) — 1/2. 
This fact was established numerically in Ref. [23j , where 
the corresponding randomness was dubbed 0(1). For 
two other types of sign disorder: random-bond Ising 
model and Cho-Fisher (CF) model [Hj, the metal- 
lic phase is either absent or emerges when W exceeds 
certain W c , correspondingly. The above two facts were 
also established numerically by studying transmission of 
finite- width (up to M = 256) stripes 0,[23|-[25|. 

In conventional quantum Hall effect, emergence of an 
isolated delocalized state at (t 2 ) = 1/2 is qualitatively 
transparent, since the Chalker-Coddington model [26| is 
a quantum version of a classical percolation, which pos- 
sesses self-duality. By contrast, there is no classical ver- 
sion of the D-class. Thus there is no qualitative expla- 
nation how a particular local sign disorder in U and Vi 
results either in metal or in insulator at large distances. 
Moreover, the fact that critical disorder in CF model is 
weakly sensitive to (M) [23|, |24[ seems counterintuitive. 
Indeed, it implies that derealization takes place no mat- 
ter how weakly the M = contours are coupled to each 
other. 



In this paper we demonstrate that above puzzles of the 
D-class derealization find a natural explanation within 
a simple real-space renormalization group (RG) descrip- 
tion. This description is developed below by adjusting 
to the D-class the approach of Ref. [27j developed for 
conventional quantum Hall effect. We focus on the CF 



model most general in the sense that its phase diagram 
in the ((£ 2 ), W)- plane contains all three (metal and two 
quantum Hall insulator) phases. 

The RG procedure [27[ prescribes how the magnitude 
and sign of effective transmission coefficient, t, evolve 
with the system size, L. Within RG we demonstrate that 
I-M transition (finite i 2 at large L) occurs at the same W c 
where the distribution of the amplitude, £, becomes sym- 
metric with respect to i = 0. For smaller W (insulating 
phase), the initial signs of U are "forgotten" with increas- 
ing L. In this sense, the I-M transition can be viewed as 
sign percolation. Remarkably, RG description appears to 
be very accurate on the quantitative level. In fact, it re- 
produces the entire phase diagram found in [23l425j. In 
addition, we find the following universal distribution of 



the conductance, G = t , of the thermal metal 



P[G) = 0.414[G(1 - G)\ 



-0.6 



(2) 



Within RG, insensitivity to (M) emerges as a result of 
high (resonant) transmission via certain local disorder 
configurations; importantly, these resonances proliferate 
to large distances. 

RG transformation. In the "clean" case, the network 
model of D-class is shown in Fig. [2j The signs of trans- 
mission and reflection at each node ensure the 7r-flux per 
plaquette. With respect to signs, scattering matrices of 
the nodes are different for two sublattices of the square 
lattice, so that S'-matrix has the form Eq. (pQ) for red 
nodes and 



(3) 



for green nodes. The RG superblock shown in Fig. [2] 
consists of four red and one green nodes. With signs in 
the S'-matrices chosen accordance with Fig. [2j scattering 
matrix of this superblock reproduces the form Eq. (pQ) 
[IK with effective transmission coefficient given by 



t 



tih(r 2 r^r A + 1) + t 2 t A {rir^r b + 1) + t 3 (M 4 + hh) 
(r 3 + rir 5 )(r 3 + r 2 r 4 ) + (t 3 + tife)(*3 + Ms) 

At the next RG step, values of U generated by Eq. (|4j), 
including the signs, should be placed into red nodes of 
renormalized lattice, which has a doubled lattice constant 

The limit W = 0. As a first test of applicability of the 
RG description we apply Eq. (jl]) to an ordered system. 
Upon setting all U equal to t and all ri equal to yl — t 2 , 
we see that t = t = l/y/2 is indeed a fixed point, corre- 
sponding to zero mass. Finite \t 2 — 1/2| plays the role of 
a finite mass, M, in a clean system. From Dirac equation 
it follows that transmission disappears at lengths exceed- 
ing tunneling length 1/M. This means that the W = 
critical exponent is v = 1. On the other hand, Eq. (J2J) 
transforms finite \t 2 - 1/2| into \i 2 - 1/2| = r\t 2 - 1/2|, 
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where 



12 



dt 2 



= 1/2 



(5) 



Elementary calculation yields c\ — c 2 — ca = £5 = y/2 — 
1, c 3 = 3 - 2>/2, so that r = 2^ - 1. Since r > 1, 
transmission disappears upon n subsequent steps, where 
n is determined by the condition, |t 2 — l/2\r n ~ 1. For 
such n the unit cell is £ = 2 n and should be identified 
with the in-gap tunneling length. From the above two 
conditions we find that f ~ |t 2 - l/2\~ v = \M\~ V , where 
v = In 2/ In r ~ 1.15, i.e., the RG value of the exponent 
is only 15% different from the exact value. 

Before introducing disorder in signs of U and n, we 
elucidate another property of t 2 = 1/2 fixed point: ini- 
tial symmetric distribution of t 2 around 1/2 not only re- 
mains symmetric but narrows upon renormalization. In- 
deed, the width of distribution after one RG step evolves 

as (St 2 ) = T?i=i c i( St i) ~ °- 7 (^;)- This means that 
the critical states at the above fixed point do not ex- 
hibit mesoscopic fluctuations, in contrast to conventional 
quantum Hall critical state. We will see below that ther- 
mal metal does exhibit strong mesoscopic fluctuations. 

Resonances. Suppose that in Eq. (j4]) all U are small. 
Corresponding Ti are close to 1. It might seem that t is 
even smaller, oc t 2 . At this point we divulge the following 
remarkable property of the transformation Eq. (|4]). Let 
us choose the following sign combination: 



ri = r 2 = r 3 = r 4 



(6) 



Substituting the above values into Eq. (j4]), we find that 
t is identically equal to ±1, i.e., resonant tunneling takes 
place. It might also seem that such resonant configura- 
tions are irrelevant due to vanishing statistical weight. 
On the contrary, we will see that these resonances play a 
key role in derealization, since they persist even when n 
are not equal. Evidence to this fact can be found in Fig. 
[3^,b. Bare values of tf are homogeneously distributed in 
the interval, 0.1 < t 2 < 0.3. For a modest randomness in 
signs, Eq. (J2J) transforms the box-like distribution into 
a three-peak distribution; the left peak corresponds to 
resonant reflection, the middle peak corresponds to ex- 
pected overall reduction of transmission, while the right 
peak comes from resonant transmission. We see that, 
with spread in the bare tf values, the portion of reso- 
nant transmission is considerable, about 8%. Figs. [3]> 
h illustrate that the competition between the resonant 
transmission peak and low-tf peaks decides whether the 
system evolves into a metal or an insulator. In the metal- 
lic phase, Figs. [3ji, [3f, and[3fr, the resonant transmission 
peak continuously broadens and grows until the distribu- 
tion becomes symmetric around t 2 = 1/2. In the insulat- 
ing phase, Figs. [3t, [3fe, and [3fe, the resonant tunneling 
peak gets gradually suppressed with the system size, and 
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FIG. 3: (Color online) RG evolution with system size, L, of 
the initial distribution of t\ (purple box) is shown for two 
disorder magnitudes, W — 0.15 (left) and W — 0.2 (right). 
Insets: (c) and (d) Distribution of the amplitude reflection 
coefficient, r^, in the domain — 1 < n < 0; (e) and (f) Dis- 
tribution of tf with ti < 0.025 removed illustrate precursors 
of insulator and metal, respectively; (g) Distribution of t\ in 
the domain, 10 -4 < tf < 10 -5 illustrates approach to the 
insulator. 



large-scale distribution flows to tf — >> 0. Details of deter- 
mination of the I-M boundary are described below. 
CF sign disorder. According to prescription of Ref. [3] 
for the isotropic version of the CF model, for a given 
positive t the values of U are chosen to be 



t, with probability (l 



-t, with probability \. 



(7) 



Correspondingly, 



y/l — t 2 , with probability (l — ^ 



w 



(8) 



-yl — t 2 , with probability 9 . 

The RG transformation for n has a form 

, _ rrr^MUk + 1) + r^(tit 2 U + 1) + r 3 (rir 4 + r 2 r 5 ) 
(r 3 + rir 5 )(r 3 + r 2 r 4 ) + (t 3 + tit 2 )(t 3 + t 4 t 5 ) ' 

(9) 

which is in accord with Eq. ((4]). With prescription Eqs. 
(O, jEj), and Eqs. (j4j), (|9j), Mathematica performs each 
subsequent RG step almost instantaneously even for a 
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W=0.05 <t 2 >=0.5 W=0.07 



FIG. 4: (Color online) RG evolution with sample size, L, of 
the symmetric distribution of 0.25 < tf < 0.75 (purple box) 
is shown for two magnitudes of disorder, W = 0.05 (left) and 
W = 0.07 (right). 



dense sample of 10 6 U and T{. We started with checking 
that for initial box-like distribution of t 2 in the interval, 
tfnin < t 2 < tmaxi the evolution depends only on the 
average, (t 2 ) — (t^ in + t 2 nax )/2, but not on the width, 
(tmax-tmin)- For particular box, t 2 min = 0.1, t 2 max = 0.3, 
results shown in Fig. [3]illustrate that distinction between 
two disorder values, W = 0.15 and W = 0.2, emerges at 
system size, L = 2 5 . At the same time, the difference 
in signs distribution of develops already at the sec- 
ond step. This is illustrated in Figs. [3t,d insets. This 
behavior is generic: memory about initial sign disorder 
in T{ is quickly forgotten in the insulating phase (all r\ 
have are positive after 3-4 steps), whereas in the metallic 
phase both signs of are completely equilibrated after 3- 
4 steps. Nevertheless, complete convergence to the sym- 
metric fixed point distribution P(t 2 ) given by Eq. ([2]) 
takes place at very large size, L = 2 25 . 

RG evolutions depicted in Fig. [3] suggest that a M-I 
transition point lies within 0.15 < W c < 0.2. By grad- 
ually shrinking this interval, one can find W c with high 
accuracy. The results for different (t 2 ), which are in gen- 
eral agreement with numerics of Refs. [23l-[25|, are shown 
in Fig. [TJ We see that the I-M boundary is approximately 
horizontal, except for the interval, 0.4 < (t 2 ) < 0.6. 
Within this interval the boundary rapidly drops to the 
value, Wtr = 0.06, which we identify with trier it ical 
point. Insets in Fig. [1] illustrate our statement concern- 
ing sign percolation. They show how the portion, ft, of 
negative revalues Eq. (j9j) evolves with the sample size. 
Equilibration of signs in metal, W > W c , implies that 
ft grows towards ft = 1/2. We see that the bigger is 
(W — W c ) the faster is the growth. By contrast, a de- 
crease of ft with L in insulator implies that all Vi become 
positive after several steps. It is seen that erasing of signs 



is more efficient at smaller W . As signs are erased, the 
magnitude, (r 2 ) approaches 1. We used the rate of this 
approach to estimate the critical exponent of the I-M 
transition. Assuming that (1 — (r 2 (L))) is a function of 
a single parameter, (W c — W) U L, we found v « 1.2 for 
initial distribution with (t 2 ) = 0.2 (red triangle mark in 
Fig. [Q). 

Tricritical point. In the RG language, trictritical point 
on the vertical axis, t 2 = 1/2, suggests that a symmetric 
initial distribution of t 2 evolves to P(t 2 ) = 5(t 2 — 1/2) 
for W < W tr , and to metallic fixed point, Eq. (|2]) for 
W > W tr - Fig. (j4]) shows the evolution of symmetric 
box-like distribution 0.25 < t 2 < 0.75 for two VK-values. 
We see that for W = 0.05 and W = 0.07 these evolu- 
tions are different. For W = 0.07 the distribution at 
L = 2 3 is flat. After a few steps the histogram bends 
down at the center and distribution flows towards metal- 
lic fixed point Eq. (j2j). For W = 0.05, the box-like 
initial distribution develops a maximum at L = 2 3 , Fig. 
2h. This narrowing suggests a flow towards the ordered 
fixed point, t = r = 1/V2. By gradually shrinking the 
VF-interval we locate the value W tr ~ 0.06, which sepa- 
rates the behavior Fig. [4J2 (maximum at the center) and 
Fig. 0]i (minimum at the center). We identify W tr with 
tricritical point 0, Q EH-El. The smaller is W the 
more pronounced is the shrinking of initial distribution 
with L. Note, however, that this shrinking eventually 
stops: at large L maximum at the center is accompanied 
by satellite peaks at t 2 = and t 2 = 1, which gradually 
take over and drive the system to metal. Similar com- 
plications were pointed out in Ref. [14|. The behavior 
of sings of ti (ri) is synchronized with distribution of t 2 
(r 2 ). Namely, as the distribution shrinks, the portion, ft, 
changes from W down monotonically. As the distribution 
turns back to metal, ft starts to grow towards ft = 1/2. 
Discussion The value, W c ~ 0.2, and the fact that it 
depends weakly on t 2 , can be inferred from the calcula- 
tion of likelihood, Vw, of resonant configurations of the 
type Eq. (|6]). These configurations occur when denom- 
inator in Eq. (|4]) is small. When all U are small and 
\ri\ close to 1, this condition requires that the product, 
(rs + rir^)(rs + ^4), is small. On the other hand, if 
both brackets are small, we will have t ~ ti, i.e., the res- 
onance is absent (due to suppression of the numerator). 
The probability that only one of the brackets is small is 
given by 

v w = 4jy(i - wf + 4w 3 (i -w)= 1 ~^ 1 ~ 2W "^ 4 . 

(10) 

It turns out that Vw is almost flat and close to 1 /2 in the 
wide interval, 0.2 < W < 0.8. This is the consequence 
of the fact that first three derivatives of Vw are all zero 
at W = 1/2. Therefore if the resonances do not prolif- 
erate at W ~ 0.2, they will not proliferate upon further 
increase of W . This explains why, at small ti, the critical 
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W c is close to 0.2. Another feature of the phase diagram 
Fig. [TJ weak dependence of W c on (t 2 ), can be traced to 
the fact that for resonant configuration Eq. © we have 
t 2 = 1 regardless of the value of r. 
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